Geometry and Spectrum of Casimir Forces 
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We present a new approach to the Helmholtz spectrum for arbitrarily shaped boundaries and 
general boundary conditions. We derive the boundary induced change of the density of states in 
terms of the free Green's function from which we obtain non-perturbative results for the Casimir 
interaction between rigid surfaces. As an example, we compute the lateral electrodynamic force 
between two corrugated surfaces over a wide parameter range. Universal behavior, fixed only by 
the largest wavelength component of the surface shape, is identified at large surface separations, 
complementing known short distance expansions which we also reproduce with high precision. 
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The famous question "Can one hear the shape of a 
drum?" was posed in 1966 by Kac, illustrating the prob- 
lem to deduce the shape of a region from the knowledge 
of its resonance spectrum . It was answered negatively 
but the difficulty to characterize the distribution of 
eigenfrequencies of the Helmholtz wave equation in arbi- 
trary geometries remains. This is particularly highlighted 
by the long-lasting efforts for chaotic (quantum) billiards 
which are described in two dimensions also by the wave 
equation 0|. The same problems occur for Casimir in- 
teractions in three dimensions p| which are of much re- 
cent interest mainly due to the advent of improved 
experimental techniques yielding the force between two 
parallel plates and between a plate and a sphere Q. 
Naturally, the question arises to what extent the Casimir 
force characterizes the shape of the interacting objects. 
A large body of theoretical work, including proximity and 
pairwise additivity approximations || , semi-classical ap- 
proaches based on Gutzwiller's formula a multi- 
ple scattering expansion perturbative methods [llj 
and recently proposed approximations from classical ray 
optics [13, has been used to compute Casimir forces. 
However, these approaches either neglect diffraction or 
are limited to small smooth deformations, and reliable 
general expressions are not known even for simple ge- 
ometries. In this Letter, we present a formula for the 
density of states which is formulated in terms of the free 
Green's function at the boundaries only. We demonstrate 
its applicability by computing the lateral Casimir force 
between two corrugated surfaces which has been studied 
in a recent experiment We find that the surface's 
shape can be deduced from the force at short distances 
whereas at large separations universality prevails. 

The electrodynamic Casimir energy of two discon- 
nected metallic surfaces S a , ct — 1,2, is determined 
by the change in the photon density of states (DoS) 
Sp(k) = p{k) — po(k) caused by moving the surfaces from 
infinity to a finite distance in vacuum where po{k) is the 
DoS for infinitely separated surfaces. Thus Sp(k) con- 
tains neither volume terms nor single surface contribu- 
tions but measures only changes in geometry by moving 



the surfaces rigidly. We consider the Helmholtz equation 



(v 2 + fc 2 )<Mx) = o 



(1) 



for a scalar field <p in the three connected regions sep- 
arated by surfaces S a on which <f> fulfills Dirichlet (for 
transversal magnetic modes, TM) or Neumann (for 
transversal electric modes, TE) boundary conditions. 
The total DoS p(k) is then given by the sum of the DoS 
for the three isolated regions. For simplicity, we have as- 
sumed that the surface geometry allows for a separation 
of the electromagnetic field into TE and TM modes which 
is possible for a large class of geometries as, e.g., for uni- 
axially deformed surfaces ■ For both types of modes 
the Casimir energy is obtained as an integral over Sp(k). 
Since the DoS is more regular along the imaginary axis, 
it is useful to shift the integration to that axis and to go 
over from a Minkowskian to a Euclidean formulation by 
a Wick rotation. Then the Casimir energy becomes 
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with Sp(q ) = -Sp(iq ). 
trace formula 



dqoqoSp(q ) 



(2) 



Our main general result is the 



Sp(qo) = -i^-Trln^M- 1 ) , (3) 

where the matrix operator A4 is given by the 
Euclidean free Green's function Go(x, x'; qo) = 
e -9o|x— x I / (4tt I x — x'|) evaluated at the surfaces 
S a only. Explicitly if the surfaces S a are repre- 
sented by 3D vectors s Q (u) with 2D local coordi- 
nates u, then M a p(u, u';q ) = G (s Q (u), s^(u'); q ) 
for Dirichlet conditions and A4 a /3(u, u'; qo) = 
c , n Q (u)9 n/3 (u') G o(sa(u),s^(u');go) for Neumann con- 
ditions with <9 nQ the surface normal derivative pointing 
into the region between the surfaces. The trace in 
Eq. © is performed over the 2D coordinates and the 
discrete surface indices. is the functional inverse 

of M. for infinite surface separation. A formally similar 
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formula for the Casimir DoS has been derived by Balian 
and Duplantier in terms of a different matrix operator 
which describes surface scatterings [T^j- In the spectral 
theory of quantum scattering an expression of the form 
of Eq.@ is known as Krein-Friedcl-Lloyd formula which, 
however, applies to the S-matrix of potential scatterers 
An important advantage of Eq.JHJ) is that it yields 
directly the regularized variation of the DoS which is 
free of distance independent divergences which would 
seriously hamper any numerical evaluations. To set an 
example for the applicability of our trace formula to the 
Casimir effect in non-trivial geometries, we compute the 
lateral Casimir force between corrugated surfaces which 
is especially sensitive to geometry. For that purpose, 
we employ a previously developed numerical algorithm 
[Til Hq which provides a fast convergent result for the 
trace. It is important to note that the lateral force does 
not arise from a change of the mean surface separation 
(yielding the normal force) but from a lateral shift of the 
boundaries, and thus requires a careful regularization of 
the energy. 

At first, we give a brief survey of the steps which 
lead to Eq. We consider the Gaussian action S = 
\ J d 4 X(V4>) 2 in 4D Euclidean space-time to quantize 
the modes of the electromagnetic field in the regions 
which are separated by the surfaces S a on which Dirich- 
let or Neumann boundary conditions hold. Path integral 
techniques with delta functions enforcing the boundary 
conditions have been used to compute the (free) energy of 
constrained systems 0, . The same approach, how- 
ever, can be also used to study correlations ^j- The 
modified correlations G(x,x';qo) = {<t>qo ( x )0-g o ( x ')) m 
the presence of boundaries can be computed exactly in 
the present case of a quadratic action. If x and x' denote 
two positions which are located both in the same region, 
one finds 

G(x,x';q ) - G (x,x';q ) = -^ du du' G (x, s a (u);q a ) 

a/3 J 

x 7W^(u,u';go)Go(x',s /3 (u');go), (4) 

where Go is Green's function in unbounded space. M~ l 
is the functional inverse, taken with respect to u, u' 
and a, f3, of the operator M. defined after Eq. J2J, 
see above. For a fixed region, the DoS on the imagi- 
nary axis is related to the Euclidean Green's function by 
p{i<lo) = (Zqo/n) J d 3 xG(x,x;q ) where the integration 
extends over the given region. We are actually interested 
in the sum of the DoS's of the three regions into which 
free space is divided by the surfaces with the bulk contri- 
bution subtracted. Since Eq. holds for every region, 
we obtain the change 5p(iqo) by integrating the r.h.s. of 
Eq. Q with x' = x over unbounded space. Explicit in- 
tegration is enabled by the simple form of Eq. J3J with x 
occurring only in the free Green's function. Performing 
the integration both with M and M.<x, and taking the 



difference of the two results, we obtain Eq. ©. 

In the following, we consider the geometry depicted in 
Fig. ma). I* consists of two surfaces with uniaxial rect- 
angular gratings along the x\ axis with equal amplitude 
a and wavelength A. The surfaces are laterally shifted 
by b and have a mean separation H, yielding a minimal 
gap 5 = H — 2a. The Casimir energy of these surfaces at 
fixed H has to be a periodic function of b, and it should be 
minimal if the surface area with minimal surface distance 
is maximal, i.e., for b = A. This leads to a lateral force 
-Flat = —d£/db. Performing the frequency integration by 
parts in Eq. J2J) with Eq. Q , we find 
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For periodic geometries as the one considered here, the 
trace can be computed with the technique introduced in 
0,0]. We Fourier transform M. a p(u, u'; qo) with re- 
spect to u and u' so that use can made of the periodicity 
along xi and translational invariance along x 2l allowing 
for the representation 



M a p (p, q; go) = (2tt) 2 5 ( P2 + q 2 ) 

OO 

x HPi+qi + ^m/X) N%f ( qi ,q 2 ;q ) , 



(6) 
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which defines the 2x2 matrices Nm that can be com- 
puted analytically and are given in |l4j for the geome- 
try of Fig^a). The trace in Eq.J^J) is more efficiently 
obtained if it is restricted to momenta q± in the inter- 
val [0, 2?r/A). This is possible after a rearrangement of 
the elements of M. so that it has block-diagonal form 
where the blocks are numbered by the continuous in- 
dex q\ in [0, 27r/A) and, for fixed qo, q 2 , have matrix 
elements B%f (qi, q 2 ; qo) = N[ 
integer indices k 7 I = — c 



k-lfal + 2irl/\,q 2 ;q ) for 
. , oo. Physically, a block 
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FIG. 1: (a) Geometry consisting of two parallel plates with 
laterally shifted uniaxial rectangular gratings, (b) Approxi- 
mate validity ranges of proximity force (PFA) and pairwise 
(PWS) approximations and sector of asymptotic universality 
for the lateral Casimir force as estimated from Figsl^landEI 
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with index q\ couples only waves whose momenta dif- 
fer from the Bloch momentum q\ by integer multiples of 
2ir/\. Thus, in analogy to Bloch's theorem, the original 
problem has separated into decoupled subproblems at a 
fixed q\ in the interval [0, 2tt/X), and the total trace in 
Eq. (J5J is given by the sum over the traces of the subprob- 
lems. This fact can be expressed by defining the function 
g(q ll q 2 ;q ) = tr(B _1 9 b i3 - S" 1 ^^) where the trace 
runs over the indices k, I, a, j3 of B, and is the analog 
of .Moo, i-C-, B for H — > oo. The lateral Casimir force is 
then given by 



■flat = o / d 1<l 



2-k/X 



dqig(qi,q2;qo), (7) 



with q = + q\. Since the matrices N m are known an- 
alytically [l4[j the same applies to B, B^, and the deriva- 
tive with respect to b can be easily computed. However, 
for a non-perturbative treatment, the inversions of B 1 
Boo have to be implemented numerically. This is enabled 
by a truncation of the matrix B at a fixed order M so 
that the function g is replaced by <?m which is defined as 
g above but with the trace running over k, I = —M . . . M 
only. F\ at follows then from a numerical integration over 
gu in Eq. for a sequence of fixed M and a subse- 
quent extrapolation to M — > oo. For the results shown 
below, we have chosen M between 13 and 37 with the 
larger M used at smaller separations S. This is physically 
consistent with the fact that with increasing separation 
smaller integer multiples of 2ir/\ around the Bloch mo- 
menta have to be considered. It should be stressed that 
the above analysis is independent of the boundary condi- 
tions which, however, change M. Thus the electromag- 
netic Casimir force is the sum of -Flat for Dirichlet and 
Neumann boundary conditions, respectively, leading to 
the results for F\ at summarized in Figs. l2l3l 

At short surface distances approximative methods can 
be employed and it is instructive to compare their pre- 
dictions to our findings. To begin with, the proximity 
force approximation |5| yields the lateral force Fpfa = 
[2£ (H)-£ (H-2a)-£o{H+2a)]/XforO < b < A/2 with 
£q{H) = — (n 2 / 720) he/ H 3 since it sums the flat plate en- 
ergy £q at the local distance normal to the surfaces. Fpfa 
changes sign at b = A/2 discontinuously. A different ap- 
proximation consists in the pairwise summation (PWS) 
of Casimir-Polder potentials. Although strictly justified 
for rarefied media only, it is often also applied to met- 
als, using the two-body potential U(r) — —{■K/24)hc/r 7 
whose amplitude is chosen as to reproduce the correct 
result for flat ideal metal plates 0. It yields the lateral 
force F PWS = J v d 3 x J v d 3 x'C/(|x-x'|) with Vi and 
V r denoting the semi-infinite regions to the left and right 
of the two surfaces in Fig^a), respectively. To compute 
-PpwS) w e have first differentiated analytically with re- 
spect to b and then performed the remaining integrals 
numerically. Fig|21 shows our results for the amplitude of 
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FIG. 2: Lateral force -Flat at b = A/4 obtained from Eq. JZJ 
for the geometry of Fig. H]a) as function of the surface gap 5 
(solid curves) . Fi at is measured in units of the normal force Fo 
between flat surfaces with a — 0. Plotted are also the proxim- 
ity force (PFA, dash-dotted curves) and pairwise summation 
(PWS, dashed curves) approximations, and the perturbative 
result F p t for sinusoidal profiles (dotted curves) 
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FIG. 3: Shape dependence of Flat on the lateral surface shift 
b at fixed distance H — 10a for different corrugation lengths. 
The dashed and the dotted curves represent the PWS and the 
full perturbative result for sinusoidal profiles with arbitrary 
H/\ [Tll |. respectively. 

the lateral force, measured relative to the normal force 
Fq between flat plates at the same H, over more than 4 
orders of magnitude for the gap 5 together with the two 
approximations. For small 6, both approximations agree 
and match our results. Beyond 6 > A/20 the PFA starts 
to fail since it does not capture the exponential decay of 
Fl a t for increasing 5. The PWS approach has a slightly 
larger validity range [cf. Fig^Jb)] and reproduces the ex- 
ponential decay. However it deviates by at least one order 
of magnitude from Fl a t for S > 2.5A, see Fig|2 Thus it 
is important that in the asymptotic limit of large surface 



4 



gaps, one can expect a universal behavior of the force 
which is independent of the detailed shape of the surface 
corrugation. Precisely this expectation is strikingly con- 
firmed when we compare our results to the perturbative 
expression for the lateral force for sinus oidally shaped 
surfaces (with amplitude ciq and wavelength A) . 

_ 8ir 3 hc a 2 f 2tt \ 27rH /\ 

which is expected to hold for a <C A <C H. By con- 
sidering the lowest harmonic of the rectangular corru- 
gation, which corresponds in Eq.JHJ to ao = Aa/ir, we 
find excellent agreement between F pt and our results for 
the geometry of Fig^a) for large distances S > A, see 
FigEl Higher harmonics describing short scale surface 
structure are irrelevant for the asymptotic Casimir in- 
teraction. This is particularly highlighted by the depen- 
dence of the force on the lateral shift b shown in Fig|3 
for H = 10a and A varying in the interval indicated by 
the bar in Fig^b). With changing A, three regimes can 
be identified. For A 3> H, the force profile resembles 
almost the rectangular shape of the surfaces, and the 
PWS approach yields consistent results. For decreasing 
A, yet larger than H, the force profile becomes asymmet- 
ric with respect to b = A/4 and more peaked, signaling 
the crossover to the universal regime for A < H where 
the force profile is sinusoidal. In the latter case, for not 
too small A/a ss 10, our results for i*| at indeed agree 
well with the perturbative result for arbitrary H/X of 
Ref.H3, cf. FigH We note that the PWS approach fails 
to predict the asymmetry of the force profile, and the 
PFA even predicts no variation at all with b in the range 
of FigEI Finally, we observe a non-monotonous change 
of the lateral force with H/X when H/a is kept fixed, see 
FigQ] After a linear increase of Fi at /F , described by 
the PFA, the force shows a maximum beyond which it 
decays exponentially. For small amplitudes, a/H = 1/40 
in Fig^J the position of the maximum at H « 0.4A is 
again in good agreement with the full perturbative re- 
sult |ll|. With increasing a/H, the maximum is shifted 
towards smaller wavelengths. 

In conclusion, we have derived a formula for the change 
of the Helmholtz spectrum by arbitrarily shaped bound- 
aries. From non-perturbative results based on this for- 
mula, we argue that the lateral Casimir force between any 
two uniaxially corrugated surfaces of equal wavelength 
and amplitude is described by Eq.JHJ) for H much larger 
than the wave length. We note that the validity range for 
this universal behavior is not fully realized in the exper- 
iment of [l3j with H/a « 10 and A/a ss 55 when a is set 
to the geometric mean of the distinct amplitudes of the 
experiment. It would be interesting to probe this univer- 
sality in experiments with different shapes at larger ratios 
H/X. We studied surface deformations with a bounded 
spectrum. Stochastic surface roughness does not have 
this feature, and we then expect a different asymptotic 
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FIG. 4: Dependence of Fiat (solid curves) on the corruga- 
tion length A at fixed mean distances H. Shown are also 
the approximations (PFA, PWS) for short distances and the 
perturbative result F p t for a <C H , X. 

behavior. We note that our approach yields also the non- 
integrated DoS and can be easily used in any space di- 
mension, and also for closed boundaries. This might be of 
importance for applications to chaotic systems as quan- 
tum billiards. At finite temperatures, the DoS has to be 
integrated with a saw teeth- like weight • Finally, ma- 
terial properties of the surfaces can be included in our 
approach in form of non-local boundary conditions |18j . 
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